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The bad-metal behavior of the iron pnictides has motivated a theoretical description in terms 
of a proximity to Mott localization. Since the parent compounds of the iron pnictides contain an 
even number of 3d-electrons per Fe, it is important to determine whether a Mott transition robustly 
exists and clarify the nature of the possible Mott insulating phases. We address these issues in 
a minimal two-orbital model and a more realistic four-orbital model for the parent iron pnictides 
using a slave-spin approach. In the two-orbital model with two electrons per Fe, we identify a 
single transition from a metal to a Mott insulator, showing that this transition must exist as a 
result of orbital degeneracy. Depending on the ratio between the inter- and intra-orbital Coulomb 
repulsions, the insulating state can be either a high-spin Mott insulator or a low-spin orbital-Mott 
insulator. In the four-orbital model with four electrons per Fe, we find a rich phase diagram for 
the metal-to-insulator transition. At strong Hund's couplings, a localization transition to a high- 
spin Mott insulator always occurs. At zero and weak Hund's couplings, on the other hand, we 
find a transition to an intermediate spin insulating state. This transition can be viewed as an 
orbitally selective metal-to-insulator transition: the transition to a Mott insulator in the xz and yz 
orbitals takes place at the same critical coupling as the transition to either a band insulator at zero 
Hund's coupling or an orbitally polarized insulator at weak but finite Hund's coupling in the xy and 
— y^ orbitals. The implications of our model studies for the physics of iron pnictides and iron 
chalcogenides are discussed. 
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I. INTRODUCTION 

The microscopic physics of the iron pnictides and re- 
lated high- Tc supcrconductorsii^ is a subject of exten- 
sive studies. The parent systems are antiferromagneti- 
cally ordered^, implying that Coulomb interactions must 
play some role. The metallic nature of these systems 
gives rise to a tendency to treat the interactions per- 
turbatively. However, various considerations have led 
to the notion that the parent iron pnictides and iron 
chalcogenides are on the verge of a Mott localization 
transition. These considerations have been based on the 
observed bad-metal properties^i^, first-principles calcula- 
tions^ii, and related analyses—"—. The bad-metal prop- 
erties are characteristic of metallic systems in proxim- 
ity to a Mott localization. The electrical resistivity at 
room temperature corresponds to a short mean free path, 
on the order of the average inter-electron spacing. The 
Drude weight seen in the room temperature optical con- 
ductivity is considerably suppressed compared to its non- 
interacting valuei^"— . Finally, relatively small changes 
of temperature induce transfers of the optical spectral 
weight extending to the eV rangei^ii^iii. The lack of 
observation^^ of any pronounced incoherent peaks in the 
high-energy electron spectrum has raised some questions 
about the incipient Mott picture, but recent microscopic 
calculationai^ have suggested that this arises from a large 
damping whose effect is enhanced by the multi-orbital 
nature of the system. 

There is some evidence for the incipient Mott picture 
from the magnetic sector as well. High-energy spin- wave- 
like excitations have been seen at the zone boundaries'^. 
In addition, the total spin spectral weight is sizable; for 



instance, it is on the order of 1 /is/Fe in CaFe2As2^. 
These properties cannot be accounted for by the elec- 
trons close to the Fermi energy alone; in particular, since 
the Fermi surfaces comprise small electron and hole pock- 
ets, the spin spectral weight coming from the electronic 
states near the Fermi surfaces will be much smaller than 
that observed experimentally. Instead, the observed spin 
excitation spectrum is more naturally associated with the 
electronic states far away from the Fermi energy, as would 
be the case in a metal close to a Mott localization. 

The incipient Mott picture arises when U jt is not too 
far away from the critical value for a Mott transition. 
(Here t refers to the characteristic bandwidth of the Fe 
3d electrons, and \J a combination of their Coulomb re- 
pulsions and Hund's couplings.) In order to further sub- 
stantiate this picture, it is important to tune the system 
into a Mott insulating state. Recently, this has been 
demonstrated^i in the iron oxychalcogenides, which con- 
tain an Fe square lattice that is expanded compared to 
the iron pnictides and iron chalcogenides. The lattice 
expansion gives rise to a narrowing of the 3d-bands and 
a concomitant enhancement of I] jt, which pushes the 
system through the Mott transition and into the Mott- 
insulating regime. 

Based on the above considerations, it is very important 
to show theoretically that a transition from a metal to 
a Mott insulator generally exists in multi-orbital mod- 
els appropriate for the parent compounds of the iron 
pnictides. This is especially so given that the number 
of 3d-electrons per Fe is even in these systems. The 
Mott transition, which has long been a subject of funda- 
mental interesli'', is studied in both one-orbital Hubbard 
model with one electron per siteS^i^, and multi-orbital 
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systems at commensurate fiUingssS. In the one-orbital 
model and paramagnetic phases, the Mott transition can 
be understood within the Brinkman-Rice picture^^, with 
an interaction-induced suppression of the coherent one- 
electron spectral weight near the Fermi energy, and the 
concomitant development of incoherent spectral weight 
away from the Fermi energy^. Focusing on the param- 
agnetic phases is advantageous for considering the Mott 
transition, since an odd number of electrons per unit cell 
means that any insulating state must be the result of in- 
teractions. The picture applies in the part of the phase 
diagram above the temperature of ordering transitions 
(typically into antiferromagnetic phases); this part of the 
phase diagram widens with the increase of magnetic frus- 
tration — . 

In our case, since the parent iron pnictides contain an 
even number of (six) 3d-electrons per Fe, an insulating 
state could in principle simply be a band insulator. The 
issue is particularly pertinent given that the Fermi sur- 
faces of the iron pnictides and related systems are small 
pockets. Indeed, the non-interacting band structure con- 
sists of non-degenerate bands (two bands are called de- 
generate only when their energies are the same at every 
point in the momentum space) and is close to that of 
a semiconductor with a small overlap in energy between 
the bottom of the conduction bands and the top of the 
valence bands; each band is far away from half-filling. It 
is a priori possible that, if an increasing U suppresses 
a metallic state, the system goes into a correlated band 
insulator first, thereby invalidating the picture of a bad 
metal on the verge of a Mott transition. In other words, 
it is a non-trivial question in the case of the iron pnic- 
tides as to whether the metallic state can be in proximity 
to a Mott transition and, if so, whether the Mott insu- 
lator resembles that for the canonical one-orbital Hub- 
bard model with one electron per unit cell. Moreover, 
given that the bandwidths are very different from band 
to band, another non-trivial question is whether a Mott 
transition, when it exists, is a one-step transition or it 
can be an orbital selective Mott transition (OSMT).— 

To address these issues, we study how a metal-insulator 
transition (MIT) in the paramagnetic phase may happen 
in multi-orbital models of iron pnictides containing an 
even number of electrons filled. Our work builds on ear- 
lier studies of multi-orbital models for transition metal 
oxides and related systems. Physics in these multi-orbital 
models is very rich. For instance, many factors, such 
as inequivalent bandwidths, crystal field splitting, and 
Hund's coupling, may affect the nature of the transition 
i^^i^. We note that there is a tendency in the litera- 
ture to work in the band representation and ignore the 
orbital characters .^^^^ In our study, we show that the 
orbital characters serve as another important factor for 
the Mott transition. 

We first study the MIT in a two-orbital Hubbard model 
with two electrons per site, with a kinetic-energy part 
given by the minimal band dispersion for the pnictides^. 
We analyze the model within a slave-spin (SS) formula- 



tio n^ 1'^^ , which has the advantage of readily capturing 
both the coherent and incoherent part of the electronic 
spectrum and, in addition, the effect of Hund's coupling. 
We show that the orbital degeneracy guarantees the exis- 
tence of a one-step metal-to-Mott-insulator transition in 
this model. Neither a band insulator nor an OSMT may 
take place. The critical coupling Uc, which is larger in 
multi-orbital systems than in the single-orbital cas o^^'^° , 
is greatly reduced by a nonzero Hund's coupling; this 
is consistent with earlier studies in other multi-orbital 
contexts4i We find that the nature of the Mott insulator 
depends on the ratio of the inter-orbital and intra-orbital 
Coulomb repulsions. It can be either a high-spin Mott 
state which is driven by the intra-orbital coupling, or 
a low-spin orbital-Mott state driven by the inter-orbital 
coupling. 

We then consider a more realistic four-orbital model 
including xz, yz, xy, and — y^ orbitals. In this model, 
the nature of the MIT depends on the strength of the 
Hund's coupling. A strong Hund's coupling stabilizes a 
high-spin Mott state on the insulator side. But when 
the Hund's coupling is either zero or weak compared to 
the crystal field splitting, we find a novel orbital selective 
MIT. The xz and yz orbitals experience a transition to a 
Mott insulator; at the same critical coupling, a transition 
to either a band insulator (at zero Hund's coupling) or 
an orbitally polarized insulator (at nonzero albeit weak 
Hund's coupling) takes place in the xy and x"^ — y^ or- 
bitals. In this case, the insulating state always has an 
intermediate spin value, even at zero Hund's coupling. 
We establish this to be a direct consequence of the dou- 
ble degeneracy of the xz and yz orbitals. 

The remainder of the paper is organized as follows. 
In Sec. II we summarize the slave spin formulation and 
introduce the two-orbital model. Sees. HI is devoted 
to the MIT in the two-orbital model at half-filling. In 
particular, we propose the concept of low-spin orbital- 
Mott state. The investigation of MIT in the four-orbital 
model is presented in Sec. IV, where a rich phase diagram 
is given. Sec. V contains some concluding remarks. 



II. MULTI-ORBITAL HUBBARD MODEL AND 
SLAVE-SPIN REPRESENTATION 

The Hamiltonian for the multi-orbital Hubbard model 
is H = Hq + Hint ■ Here, Hq is a non- interacting tight- 
binding Hamiltonian with the general form 

af^u i,(T i.acr 

in real space, where c?;^,^. creates an electron on site i, in 
orbital a, and with spin a. Sa is the on-site potential of 
orbital a that incorporates the crystal field splitting, /i is 
the chemical potential determined by the electron filling. 
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The interaction reads 
Htnt = U ^nia^riiai + [U' - J/2) ^ riianip 

ia i,o</3 
- J ^ [2Sia • Si^ - (rfLt^aiO'wrfw + h.C. 

(2) 

in which U{U') denotes the intra- (inter-)orbital Coulomb 
repulsion and J the Hund's coupling. These three param- 
eters satisfy U' = U — 2J based on the consideration of 
rotational symmetry. It is strictly satisfied for an isolated 
atom in free space where all the d-orbitals are degener- 
ate and the Coulomb potential has the full rotational 
symmetry, and is assumed to be also valid in solids4^ 
We will adopt this widely used relation unless otherwise 
specified (see discussion in Sec.IIIB). The spin operators 

are Si„ = \ Yjaa' dl^a'^Ta' dima' , where r = {t"" ^t^ ,t'') 
are the Pauli matrices. 

It is now generally accepted that the degenerate d^z 
and dyz orbitals contribute most to the low-energy 
physics of the parent iron pnictides. Hence in this pa- 
per, we will first consider a two-orbital model intro- 
duced in Ref. [s^ The simplicity of the model makes 
it easier to bring out some essential insights, which will 
also be instructive for the understanding of more real- 
istic models with a larger number of orbitals. Defin- 
ing -^L = {d\^xa,d^y„), we have Hq = Y.\,a ^lAi'^ + - 
-I- £_r^ -I- SxyT^]'ip-kcr, where £+, e_, and £xy are the 
intra- and inter-orbital hopping matrices in the momen- 
tum space. In the notation of Eq. [U the orbital degen- 
eracy requires an a-indcpendent which can be set to 
zero. We also notice that the tight-binding Hamiltonian 
in this model is symmetric under the orbital interchange 
xz ^ yz. The parent compound has a half-filling, i.e., 
two electrons per site. 

We study the MIT in this two-orbital model us- 
ing the SS formulatiouj^i^^ This formulation involves 
a much smaller number of slave fields compared to 
the atomic-configuration-bascd slave-boson representa- 
tion of Ref. H^, and more readily treats the full Hund's 
coupling compared to the slave-rotor representation of 
Refs. m and 13. Here, a slave quantum S = 1/2 
spin is introduced on each site for each orbital and spin 
degree of freedom: di^a 2Sf^^fiaa- The Hilbert 
space spanned by the SS and auxiliary fermions arc lim- 
ited to the physical part by imposing the constraint 
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niaa on each site. The SS formulation 



handles the electron interactions by rewriting Hint in 
terms of SS operators. The density-density interactions 
in Eq. [5] (including the Ising-type Hund's coupling) are 
easily handled by the z-component of the SS operator. 
The spin-flip part of the Hund's coupling and the pair- 
hopping term arc approximately treated by substituting 
the fermion operators by SS operators that have the same 
effect on the SS quantum numbers of the Hilbert space, 

viz. —Jj2i '5'i^|^t'^iU'^i^4.'^i2T ~ '^ilt'^ili'^i2^'^i2t ' 

This approximation should capture the qualitative 



physics because the SU(2) spin-rotational symmetry of 
the Hund's coupling is still preserved in the Hilbert space 
spanned by the slave spins 4^ 

The SS formulation is treated at the mean-field (MF) 
level by fully decoupling the SS and auxiliary fermion 
operators via a saddle-point approximation. This leads 
to two decoupled MF Hamiltonians for the SS and the 
auxiliary fermions: 



+ h^SU + 1/2) + Hi 



(3) 



a.(3<7 iu 



(4) 




U -U' 



- J 



H.c. 



q+ c- q+ q- 



(5) 



and haa is a Lagrangian multiplier taking account for 
the constraint. To solve these two Hamiltonians, we 
further apply the mean-field decomposition to the term 

S^aaSi+u.lSa , i.e., ■S'i^Qcr'S'f+iy.^cr ~ (Staa) ^i+u.fia + 

^i:a.(^i+.,^^a> - {Sia.Sl,,.p.), ^ud assumc to be 

site independent. The Hamiltonian for the SS operators 
is then reduced to (up to a constant) 



H!iF = E [K^^-^laa + haASU + 1/2)] + H, 



s 

int ' 



(6) 

where = »EB{S^pJEiMpfLji+-^^P-)- The 

quasiparticles near the Fermi level are described by the 
auxiliary fermion Hamiltonian H-^ . Introducing the 
quasiparticle spectral weight = 4:{S^^)'^, is writ- 
ten as 

which has a similar form as Hq, with the hopping tj^^ 

renormalized to ^ ZaZ^t"^^. In practice, Eq.|6]and Eq.[7] 
are self-consistently solved by iteratively determining the 
parameters haa and K^c; (hence Z^). The metallic be- 
havior corresponds to the Bose condensation of the slave 
spins, which is marked by a non-zero Z^,- The Mott insu- 
lating behavior of orbital a is then identified by Z^ — 0. 
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FIG. 1. (Color online) Evolution of the spectral weight Z 
in the two-orbital model with J = 0, f/' = f/ in (a) and 
J/U = 0.2, U' /U = 0.6 in (b). D = 12, is the full bandwidth 
of the non-interacting band structure. 
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FIG. 2. (Color online) DOS in the two-orbital model at vari- 
ous U values with the same model parameters as in Fig. [T] 



III. RESULTS FOR THE TWO-ORBITAL 
MODEL 

A. Mott transition in the two-orbital model 

In the two-orbital model for iron pnictides, the non- 
interacting tight-binding Hamiltonian contains an inter- 
orbital hopping Exy , and cannot be diagonalized by a uni- 
form (k-independent) orbital rotation. This leads to two 
non-degenerate bands which has asymmetric local den- 
sity of states (LDOS) with respect to the chemical poten- 
tial. The Fermi surface consists small electron and hole 
pockets, indicating that both bands are far away from 
half-filling. It is then important to ask whether, when 



the interactions are turned on, the system undergoes a 
transition to a Mott insulator, or the band degeneracy 
can be fully lifted so that a correlated band insulator 
may be stabilized. If a Mott transition does take place, 
it is also interesting to ask whether the transition to the 
Mott insulator is through a single MIT or via an OSMT. 

Most of previous studies on the MIT in multi-orbital 
systems work in the band presentation, assuming a di- 
agonalized band structure with zero inter-band hopping. 
This oversimplified treatment neglects the orbital char- 
acter of the model. Within this representation, the only 
way to generate two bands with asymmetric LDOS is to 
introduce a finite crystal field splitting which breaks the 
orbital degeneracy in the model. In this study we choose 
to work in the orbital representation, which takes into 
account the full orbital characters of the model. In this 
orbital representation, the symmetry of Hq and Hint un- 
der interchanging the xz and yz orbitals guarantees that 
both orbitals contribute equally to the band structure. 

The above considerations lead to two important con- 
sequences. First, rixz.a- = nyz.a- = 1/2 for any value 
of U, i.e., each orbital is exactly at half- filling. So a 
Mott transition is possible at finite U. Second, Z^z = 
Zyz ~ Z, and an OSMT cannot take place even though 
the widths of the two bands are very different. At finite 
U, the hopping parameter fj^^^ is uniformly renormalized 
to y/ZiZ^ti^ = Zt\^^ since the crystal field splitting is 
zero in this model. Hence the dispersion ek (the Fourier 
transformation of is also normalized to Ze^. We 
thus see that the topology of the band structure is un- 
changed by the interactions. At T = 0, the chemical 
potential is determined by /0(/i — Z(.-\^ = 2iVs, where 
0(a;) is the Heaviside step function, and TVs the number 
of sites. Taking Z = 1 and /i = at [/ = 0, we see that 
/ 0(//-Zek) = / 6(Z(/io-ek)), implying that /i = Z/iQ. 
The Fermi surface is determined by Z(^o ~ Ck) = 0. 
Therefore, the Fermi surface of the interacting system is 
identical to the one in the non-interacting system. Fur- 
thermore, the filling factor of each quasi-particlc band 
does not depend on the interaction and must be identi- 
cal to the value in the non-interacting case. Therefore, 
in the presence of nonzero interaction, the quasi-particle 
bands are still partially occupied, and a band insulator 
never emerges. The system stays in the metallic state un- 
til the quasi-particle spectral weight Z = 0, where both 
bands go through a transition to a Mott insulator. Our 
argument generally applies to any system with degener- 
ate orbitals. 

The above analysis is supported by the full solution to 
the slave-spin mean-field (SSMF) equations. Fig.[T]shows 
the evolution of spectral weight Z with increasing \J . For 
both J = and nonzero J, Z drops to zero at finite J7, in- 
dicating a Mott transition. At J = 0, this takes place at 
UciD = 2.66, where £) = 12 is the full bandwidth of the 
two-orbital model (in the unit of t\ , the nearest-neighbor 
intra-a;z-orbital hopping along the x direction). The crit- 
ical coupling is reduced to Uc/D = 1.49 at J/U = 0.2. 
As a complementary method, we have also applied the 
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slave-rotor formulations^ to the two-orbital model with 
J=0 and the same values for the other model parameters. 
The results (not shown) are qualitatively the same to the 
SSMF results at J = 0: a Mott transition takes place at 
Uc/D « 2.0. 

The Mott transition is best seen in the variation 
of the spectral function with increasing U. We cal- 
culate the spectral function by convoluting the SS 
and auxiliary fermion Green's functions: G^'^(k, w) = 

*/d..'{G>|(c.')G:^^(k,a.-a.')+G^^^(a.')G</(k,a;-..')}, 
where is the retarded Green's function, G^p{t) = 

-i{S^{t)S-^m, and G</(i) = tifjMUt)) and a, (3 
denote orbital indices. Note that in the above expres- 
sion the SS Green's function is independent of k. This 
is a consequence of the MF approximation. Expressing 
the SS Green's function using the Lehmann representa- 
tion and taking into account that the Hamiltonian of /- 
fermions describes free fermions, the spectral function is 
written as 

aX(T n,m 

x{e~^^-{l-ni^) + e-P^'^nQ, (8) 

where i?„ and \n) are eigenenergy and eigenvector of SS 
Hamiltonian, Z = J2n^~'^^"' ^rim = En - E,n, and 
n{ = l/{e^^^ + 1). The matrix Ak diagonalize the /- 
fermion Hamiltonian Hjjp with eigenenergy eAk. Taking 
m = n = (|0) denotes the ground state of Hfjp) in 
Eq.[51 one sees that the coherent part of the spectral func- 
tion is normalized by a factor of Z since Z = |(0|S'2£^|0)p. 
As mentioned above, an advantage of the SSMF is that 
the incoherent part is also accessible. At low temper- 
atures this comes from terms with n = 0,m 7^ 0. The 
LDOS is calculated from Eq.|8l As shown in Fig.[2l when 
U is increased from zero, one sees clearly that the coher- 
ent part is renormalized by Z . There is a significant spec- 
tral weight transfer to the incoherent part in the metallic 
phase. Our results provide a natural explanation of both 
the renormalization of the coherent bands and the ap- 
pearance of the incoherent spectral weightsi^i^ within a 
unified framework. At J7 > Uc, the coherent peak van- 
ishes, signaling the Mott insulator state; the incoherent 
parts, at the same time, develop into the lower and upper 
Hubbard bands. 

As discussed in previous studiesj^ a finite Hund's cou- 
pling may strongly affect the Mott transition. We show 
this effect in the two-orbital model by presenting the J- 
U phase diagram in Fig. [31 It is seen in Fig. [3] that 
Uc is rapidly reduced with increasing J/J7; this is also 
illustrated in Fig. [TJ The reduction of Uc can be un- 
derstood by solving the mean-field Hamiltonian in Eq. [5] 
at an infinitesimal K {K = Kaa)- We first diagonalize 
Hint and the term including S?^^ in Eq. [6l and label the 
eigenstates as \ni')^ where n is the electron occupation 
number, and v denotes the degenerate multiplets. We 
then treat the term KH' = K X^ao- ^aa pcrturbatively. 



To the first-order va K, Uc can be obtained by solving 
1/e = 2£. Here e = l/Ns^te Y^ap {flacrf^aa) is the av- 
erage kinetic energy for the noninteracting system, and 
£ is the lowest eigenvalue diagonalizing the matrix M, 
where M^, = Y.n^2,x{MH'\nX){n\\H'\2u) / [E2 - En) 
and En is the eigenenergy of state \nv). M^i, is non- 
zero only when n = 1 or n = 3. For either n value, 
E2 ~ En = -A2/2 where A2 = J7 + J is the Mott gap 
of the two-orbital model at half-filling. Hence Uc is de- 
termined by A2 cx |e|, i.e., Uc cx |e|/(l + J/U). This 
clearly indicates that Uc decreases with increasing J/U 
ratio. Similar behavior has also been discussed in a three- 
orbital modeli^ 

Interestingly, we find that Uc is reduced more sig- 
nificantly for the Ising-type Hund's coupling (i.e., in 
the absence of the spin-flip and pair-hopping terms in 
Eq. at large J/U ratio, as shown in Fig. [31 This is 
quite consistent with the results in previous studies^ 
and can also be understood from the above perturbation 
theory. When J is small, for either Ising-type or full 
Hund's coupling, all six configurations associated with 
two electrons occupying the two orbitals, denoted by |2/i) 
with ^ ranging from 1 through 6, are nearly degener- 
ate, and are strongly mixed in the (perturbed) ground 
state. But when the full Hund's coupling J is large, 
only the triplet configurations in \2fi), shown in Fig. [31 
contribute most to the ground state. For the Ising-type 
Hund's coupling, the ground state only strongly mixes 
the doublet: | t)xz\ t)yz and | Dy^. More config- 

urations mixed in the ground state correspond to more 
scattering processes between the nearly degenerate |2/i) 
and \2v) states, which promote a larger kinetic-energy 
gain in the metallic phase, thereby favoring the metal 
over the Mott insulator. Following the perturbation the- 
ory, Uc = 12|e|/(l + J/U) « 2.7i:»/(l + J/U) for in- 
finitesimal J/U; here, a nearly degenerate perturbation 
is used involving all the six low-energy multiplets. At 
sufficiently large J, Uc = 8|e|/(l + J/U) for full Hund's 
coupling and Uc — 4|e|/(l -I- J/U) for Ising-type Hund's 
coupling; here, the degenerate perturbation respectively 
involves three and two lowest multiplets for the two 
cases. These expressions are qualitatively — consistent 
with the numerical results in Fig. [31 At J/U < 0.01, 
Uc « 2.6615/(1+ J/C/) for both full and Ising-type Hund's 
couphngs; While for J/U > 0.03 the largest reduction of 
Uc is found in the Ising-type Hund's coupling. It is inter- 
esting to note that the effect of Hund's coupling on the 
value of Uc is quite similar to the effect of having more 
orbitals in degenerate multi-orbital Hubbard models4^ 
Indeed the underlying physics is related: in both cases, 
Uc is higher when the ground state mixes more nearly- 
degenerate configurations; as already mentioned, involv- 
ing more configurations helps stabilizing a metallic state 
by lowering the kinetic energy. 

Another difference from the J = case is that the MIT 
becomes discontinuous at nonzero J. It is especially sig- 
nificant for the Ising-type Hund's interaction. For the 
full Hund's coupling, the discontinuity of Z is only signif- 
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FIG. 3. (Color online) Phase diagram in the J—U plane of the 
two-orbital model at half-filling. Here we have taken U' = U— 
2J . The blue circles (red diamonds) give the phase boundary 
of the metal-to-Mott-insulator transition for the case of full 
(Ising-type) Hund's coupling. 



leant at small J values. For J /U > 0.2. the discontinuity 
is rather small (see Fig. [1}, and it is hard to distinguish 
the transition from a continuous one. These results from 
our SSMF calculation are consistent with the DMFT re- 
sults on the effect of finite Hund's couplings in multi-band 
Hubbard models with degenerate bands4i 



B. High-spin-Mott vs. low-spin orbital-Mott state 



High- Spin Mott: 

|^^>+|-i-^> U'-J 

Low-Spin Orbital-Mott: 

1^— >-|— ^> U-J 



FIG. 4. (Color online) Illustration of the high-spin Mott and 
low-spin orbital-Mott states in the two-orbital model at the 
atomic limit. In this model, the high-spin Mott state has S — 
1 and the ground-state energy U' — J; the low-spin orbital- 
Mott state has S = and the ground-state energy U — J. 

As seen in Fig. [51 the Mott transition at half-filling in 
the two-orbital model is quite similar to the Brinkman- 



Rice picture of the one-orbital case. It is then interesting 
to see whether the Mott insulating state is similar to that 
of the one-orbital case. For J > 0, the Mott insulator of 
the two-orbital model is an S" = 1 inter-orbital triplet 
state. This state is characterized by a vanishing on-site 
double occupancy in each orbital, and is the two- 

orbital analogue of the Mott insulator in the one-orbital 
model. We denote the triplet state as the high-spin Mott 
state in the two-orbital model. 

The Mott insulating state at J = in the two-orbital 
model is somewhat different. It mixes the high-spin 
(triplet) states with low-spin (singlet) configurations, 
which arc degenerate when J = 0. Therefore, the insu- 
lating state has a finite double occupancy (sec Fig. 
This implies that a spin-singlet state might be stabi- 
lized in some parameter regime. Note that the relation 
U' ~ U — 2J puts a strong constraint on the parame- 
ters, which may limit the ground state configurations to 
a relatively smaller subset. To fully study all the possi- 
ble ground state configurations, in this section, we relax 
the above constraint so that U' becomes a free parame- 
ter independent of U and Indeed, by studying the 
Hamiltonian in the atomic limit, one sees that for any 
J > the inter-orbital triplet state is only stabilized at 
U' < U. When U' > U, the ground state is an orbital 
anti-symmetric spin-singlet state 1/a/2(| ti; 0) ~ |0;ti)): 
as shown in Fig.|4j^ Though this state shares some char- 
acters as a band insulator, such as finite double occu- 
pancy and spin singlet, it is still a Mott insulator because 
the orbital degeneracy is preserved by the Hamiltonian. 
This can be immediately seen by noticing that each or- 
bital is at half-filling and the spectral weight is zero in 
this state. Note that U' > U is not enough to drive the 
system to a band insulator because the existence of the 
Mott insulator is guaranteed by the orbital degeneracy. 
To distinguish this Mott state aX U' > U from the high- 
spin (triplet) Mott state, we will denote it as low-spin 
orbital-Mott state. 

To see the difference between the high-spin Mott and 
low-spin orbital-Mott insulators, we study the MIT in 
the two-orbital model at several different U' /U ratios and 
show the results in Fig. [S] We find the MIT is discon- 
tinuous for general U, U\ and J values except for J = 
and U = U' , where a continuous transition is observed. 
Moreover, the behaviors in the metallic state at differ- 
ent U' /U ratios are quite similar. This is not surprising 
since the electron hopping mixes all configurations. Be- 
sides the double occupancy, the difference between the 
high-spin Mott and low-spin orbital-Mott states can also 
be seen by both the average value of the total spin oper- 
ator = {J2a ^-nd the orbital correlation function 
Cy = {{ni - l){n2 - 1)). (We have defined the orbital 
indices 1 = xz, 2 = yz.) In the high-spin Mott state, 
(S^) = 2 and C^'^ = 0. By contrast, in the low-spin 
orbital-Mott state, (S^) = and = -1. Ah these 
are consistent with the numerical results shown in Fig. [S] 
At U' = U, the ground state is a mixture of the two Mott 
states, hence both C^'^ and (S^) take intermediate val- 
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FIG. 5. (Color online) (a): The evolution of Z at half-filling 
with J/U = 0.2 and different U' /U ratios showing the transi- 
tion to different Mott states. (b)-(d): The evolution of (n^n^), 
Co'^, and (see text for the definitions of these quantities) 
for the same set of parameters. 



ues. 

For a fixed J/U ratio, wc find the critical coupling Uc 
for the MIT is the largest for U = U' . It is slightly re- 
duced for U' < U but greatly decreased when U' > U. 
We show that this non-monotonic behavior of Uc is re- 
lated to the different nature of the Mott states. Note that 
the two Mott states a±U' <U and U' > U have different 
Mott gaps. In the high-spin Mott state, A2 = U +J, is in- 
dependent of U'. On the other hand, the Mott gap in the 
low-spin orbital-Mott state is A2 = 2U' — U + J, which 
increases with U'. Therefore, Uc decreases drastically 
when U' /U increases from 1, but is almost insensitive to 
U' /U for U' < U. These considerations allow us to un- 
derstand the numerical results given in Fig. [S]^ The fact 
that Uc ai U' <U is smaller than Uc at U' = U can be 
further understood by the different ground-state degen- 
eracy in the two states. From Sec. IIIA we know that a 
higher ground-state degeneracy increases the effective ki- 
netic energy gain, thereby enhancing the stability of the 
metallic state. The ground state is three-fold degenerate 
when U' < U but is four-fold degenerate when U' = U 
(six- fold if further J = 0). Thus, Uc is the largest for 
U = U', and is slightly reduced when U' < U. 



The orbital-Mott state exists only in systems with de- 
generate orbitals and U' > U. When the orbital degen- 
eracy is broken by a crystal field splitting, it is unstable 
toward either a band insulator, or more generally as will 
be discussed in the next section, an orbitally polarized 
insulator. 
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xy 
xz, yz 



2 2 
x-y 



FIG. 6. (Color online) Three candidate ground states in the 
atomic limit in the four-orbital model. The S = 2 state is an 
analogue of the high-spin Mott state in the two-orbital model. 



IV. RESULTS FOR THE FOUR-ORBITAL 
MODEL 

We have so far shown how the Mott transition takes 
place in the minimal two-orbital model for parent iron 
pnictides. It has been recognized that, to more realisti- 
cally reproduce the electron structure of the iron pnic- 
tides, all the five Fe orbitals need to be included in the 
tight-binding dispersion^i^. This raises the question as 
to whether our main results for the two-orbital model 
are applicable to the more realistic models with a larger 
number of orbitals. In the five-orbital model of Ref. ISSl . 
the electron filling is about 0.8 per site per spin in the 



3z" 



orbital, but very close to 0.5 per site per spin 



in all other four orbitals. The Sz^ 



orbital hardly 



contributes to the band structure near the Fermi level. 
These suggest that one may study a model including only 



xz, yz, xy. 
3z2 



and 



y orbitals by assuming that the 



orbital lies far below the Fermi level and is fully 
occupied. Taking model parameters of Ref. [H but keep- 
ing only those four orbitals gives rise to Fermi surfaces 
that almost identical to the ones for the five-orbital cases. 
Hence the four-orbital model represents a good approx- 
imation to the five-orbital one. For simplicity, here we 
study the MIT in this four-orbital model. 

Wc argue that the main results for the two-orbital 
model still hold in the four-orbital model. Though the 
orbital degeneracy of other orbitals are lifted, the xz and 
yz orbitals are still degenerate. In the atomic limit, with 
four electrons occupied, the ground state may be either 
a high-spin S = 2 state, or an intermediate-spin S' = 1 
state, or a low-spin S' = state, as shown in Fig. |6l In 
either case the xz and yz orbitals are half filled just as 
in the two-orbital model. Therefore, a Mott transition 
similar to the two-orbital case is expected. 
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FIG. 7. (Color online) SSMF results for the four-orbital model 
at J = 0, showing evolution of Z (in (a)), average electron 
filling per site per spin (in (b)), average of total spin (in (c)), 
and inter-orbital correlation between xz and yz orbitals (in 
(d)). 
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FIG. 8. (Color online) SSMF results for the four-orbital model 
at 3 IV = 0.007. 



In Fig. [7] to Fig. [TT] wc show the results from SSMF 
calculation, in which the full Hund's coupling and the 
constraint [/' = [/ — 2 J are taken into account. For both 
J = and J > 0, an MIT is observed. For J = 0, 
Uc = 12.5eV. A nonzero J may significantly reduce Uc- 
At J/U = 1/4, Uc is reduced to 3.82el^. Interestingly, we 
find that the nature of the insulating state, and hence the 
MIT. is significantly affected by the competition between 
the Hund's coupling J and the crystal field splitting A 
between the xy and — orbitals. For J/U > 0.009 
the insulating state is a high-spin state with 5 = 2 as 
illustrated in Fig. [51 This state is analogous to the high- 
spin Mott state discussed in the two-orbital model. One 
may directly check from Fig.|n{d) and Fig.lTUTc) that the 
inter-orbital correlations are very small in the metallic 
state, and vanish in the insulating state, with a behavior 




12 3 4 

U(eV) 



FIG. 9. (Color online) SSMF results for the four-orbital model 
at J/U = 1/4. 



similar to the two-orbital case at finite J. 

The transitions at small and zero J/U ratios are of 
special interest. As shown in Fig. [Zl^c), the insulating 
state at J = has an intermediate spin between S = 
and S = 1. We call this as ISl state, which is shown 
in the phase diagram, c.f. Fig. [TT] This transition to 
the ISl state can be understood as follow: at J = 0, 
the xy orbital is empty and the x^ — y^ orbital is fully 
occupied in the insulating state due to the crystal field 
splitting. Hence ai U ^Uc these two orbitals undergo a 
transition to a band insulator. On the other hand, the 
degenerate xz and yz orbitals are at half-filling, and a 
Mott transition takes place at the same Uc value. Since 
at J = all the six n — 2 configurations in xz and yz 
orbitals are degenerate, the Mott insulator is a mixture 
of 5* = and 5 = 1 states. This gives the ISl state 
an intermediate spin value. The transitions to the band 
insulator and Mott insulator are reflected in the behavior 
of the intcr-orbital correlation functions Cq^ and Cq"*^, 

where Cq'*^ = ((7i„ - l)(n^i - 1)) fora,/3 = 1,2,3,4. Here 
1,2,3,4 denote xz, yz, x^ — y^ and xy, respectively. At 
J = 0, C^q" jumps to — 1 at C/c, indicating a transition to 
a band insulator. But Cq^ > — 1 even in the insulating 
phase, signaling that the insulating state is a mixture 
of S* = 1 and 5 = Mott states, just as in the two- 
orbital model at [/ = [/'. In general, the Mott transition 
in the xz and yz orbitals may take place at a U value 
different than Uc- But in this four-orbital model, Uc for 
the transition to the band insulator in the xy and x^ — 
y^ orbitals is larger than the critical value for the Mott 
transition in two degenerate orbitals. Hence when the 
other two orbitals become band insulators at Uc and thus 
are decoupled from the xz and yz orbitals, the xz and yz 
orbitals enter the Mott insulating state immediately via 
a first-order transition. Wc call this special transition at 
Uc an orbitally selective MIT. This transition is different 
from the OSMT in that it takes place at a single Uc, with 
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different orbitals entering different insulating states. 
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FIG. 10. (Color online) The evolution of inter-orbital cor- 
relation between the xy and — orbitals, Cq*, in the 
four-orbital model at J = 0, J/U = 0.007, and J/U = 1/4, 
respectively. 



We find that the orbitally selective MIT extends to 
nonzero J values up to J/U ~ 0.009. For nonzero J 
in this range, the state on the insulator side has an in- 
termediate spin value S ~ \. We denote this state as 
the IS2 state. Shown in Fig. |S] is the transition to the 
IS2 state in the case of J/U = 0.007, for which the IS2 
state is found to be stabilized for 12eV U < 15.2eV. 
Similar to the ISl state at J = 0, the IS2 state is sta- 
bilized because different orbitals undergo transitions to 
different insulating states: an 5 = 1 Mott insulator for 
the xz and yz orbitals, and an S* = orbitally polarized 
insulator for the xy and — y^ orbitals. The orbitally 
polarized insulator can be understood by studying the 
atomic limit of an effective two-orbital model including 
xy and x^ — y^ orbitals.— In this effective model, the 
pair-hopping term couples the two states | '\'i)xy\0)x^-y^ 
and \0)xy \ t4-)a;2-y2. When A > 2-\/2J, the ground state 
is a spin singlet cos0|O)^y|J\l,)^2_j^2 -sin 6*1 ti)xy\0)x^-y^ , 
where tan 6* = a/(A/J)2 + 1~A/J. One may check that 
n-xy ~ sin^ 0, and 713.2 _j^2 ~ cos^d. In general, each or- 
bital is only partially occupied, and n^y ^ nx2^y2, which 
is clearly shown in Fig. [Sljb) within 12eV < U < 15.2eV. 
Hence this state is different from either a band insulator 
or a Mott insulator, and is denoted as orbitally polarized 



insulating state. From Fig. [TOl we see that the inter- 
orbital correlation in the orbitally polarized state behaves 
as same as in the band insulator. To understand this, 
note that there are two special limits in the orbitally po- 
larized insulator. First, when J = and A > 0, sin 61 = 
and this state describes the band insulator with fully oc- 
cupied x'^ — y^ orbital and empty xy orbital. The other 
limit appears when A = but J > 0. This leads to 
sin 9 = cos 6 = 1/ \/2, and the state is identical to the 
low-spin orbital-Mott state discussed in Sec. IIIB. Thus 
the band insulator and the low-spin orbital-Mott insula- 
tor are adiabatically connected by the orbitally polarized 
insulating state. 

Further increasing J in the effective two-orbital model 
for xy and — y^ orbitals eventually leads to a low- 
spin to high-spin transitioui^ The ground state manifold 
changes from the 5* = orbitally polarized state to the 
5 = 1 Mott state. In the four-orbital model, this corre- 
sponds to a first-order transition in the insulating states 
from the IS2 state to the high-spin Mott insulator, which 
takes place at A = 2\f%J . Since in this study we assume 
J is proportional to J7, the low-spin to high-spin transi- 
tion is accessible by increasing U . In Fig. [5] we identify 
that it takes place at ?7 « 15.2eF for J/U = 0.007: the 
total spin jumps from S" = 1 to 5 = 2, and the filling 
factors of xy and x^ — y^ orbitals rapidly converge to 
half-filling. 

In Fig. [11] we show the J-U phase diagram for the 
four-orbital model. For J/U > 0.009 the MIT takes place 
between a paramagnetic metal and S = 2 high-spin Mott 
insulator. In this regime, the phase diagram is similar to 
the one in the two-orbital model shown in Fig. [3] The 
main difference from the phase diagram in Fig. |3] lies 
at U ^ A and J ^ A. In this regime we find two 
intermediate-spin states: the ISl state at J = and 
the IS2 state at finite J. The boundary between the 
IS2 and the 5 = 2 Mott insulator, J/U = A/2^/6[/, is 
determined by solving the Hamiltonian of the four-orbital 
model in the atomic limit. Given that A = 0.52eV in this 
model, and assuming that Uc for the MIT at J/U <C 1 
stays the same value as J = 0, we may estimate the 
tricritical point where the transition between the IS2 and 
the high-spin Mott state and the MIT meets. It is located 
at J/U w 0.0085; which is consistent wit h the numerical 
result J/U « 0.009. 

It is very important to check how the Fermi surface 
in the metallic state evolves in the presence of interac- 
tions. In the noninteracting case, the Fermi surface of 
the four-orbital model consists of two hole pockets cen- 
tered at (0, 0) point in the one-iron Brillouin zone and 
two electron pockets centered at (7r,0) and (0,7r) points, 
respectively. We have checked that for a large portion of 
the metallic regime in the phase diagram, i.e., to the left 
of the dotted line at J7 « 10 eV in Fig. [TTJ the Fermi 
surface of the interacting system has the same topology 
as that for the non-interacting system. In this regime, 
even near the high-spin Mott insulator, there are only 
very tiny changes in the size and shape of the hole and 
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FIG. 11. (Color online) Phase diagram in the J-U plane of 
the four-orbital model. The thick solid (black) line shows the 
phase boundary between the metallic and insulating states. 
The blue dashed line at J = refers to the intermediate-spin 
insulating state ISl, and the thin solid (red) line refers to the 
low-spin to high-spin transition between the IS2 state and the 
high-spin Mott state. The dotted line at ?7 « 10 eV separates 
the metallic phase into two regimes. To the left of this line, 
the Fermi surface topology is identical to the non-interacting 
case, while in the regime between this dotted line and the 
thick solid line, the topology of the Fermi surface changes 
(see text). Inset is a closer view of the lower- right corner of 
the phase diagram. 



electron pockets. This is similar to the result of the two- 
orbital model, in which the Fermi surface is always iden- 
tical to the one in the non-interacting system. But in the 
four-orbital model, we find that the topology of the Fermi 
surface may change when the system is close enough to 
the intermediate-spin states (in the regime between the 
dotted line and the thicker solid line in Fig. [TT|). In this 
regime, as U is increased, two additional small electron 
pockets centered at (tt, 0) and (0,7r) points and one ad- 
ditional hole pocket centered at (tt, tt) point may appear. 
Such a change in the Fermi surface topology primarily 
reflects the difference in the electron occupancies of the 
different orbitals compared to the non-interacting case 
(cf. Figs.[7t),|8b,|9lD), which is a precursor to the orbitally- 
selective MIT occurring in the small J /U regime. 

MIT has also been discussed in spin density wave 
(SDW) calculations in a number of multi-orbital models 
for iron pnictidesi^"— An intermediate-spin insulating 
state is recently found in a five-orbital model within the 
SDW MF theory,^ Though a full study of the MIT in the 
magnetically ordered state using SSMF theory is beyond 
the scope of this paper, we find it quite interesting to 
compare the intermediate-spin states found in our study 
with those in Ref. [H in the atomic limit. In both works, 
the intermediate-spins states are found when U /t is large 
but J /U is small, indicating that these states all origi- 
nate from the interplay of the crystal field splitting and 
Hund's coupling. But there arc some differences. First, 



in Ref. [H, a S' = state violating the Hund's rule is 
stabilized at J = 0. This state originates from the lifting 
of orbital degeneracy between the xz and yz orbitals by 
magnetic ordering. But In our case, the orbital degen- 
eracy between the xz and yz orbitals stabilizes the ISl 
state with a non-zero spin value at J = 0. A second 
difference lies in the way treating the electron-electron 
interaction term Hint '■ we consider the full interaction in 
our SSMF calculation; while in Ref. [H, the spin-flip and 
the pair-hopping terms are neglected from the MF ap- 
proximation. This results in different insulating states: 
the IS2 state in our paper consists an orbitally polarized 
insulator in the xy and x"^ — y^ orbitals, but the bands 
with xy and x^ — y^ orbital characters in the 5'^ = 1 state 
in Ref. [s^ are simple band insulators. 

An interesting question is whether an OSMT may take 
place when more than two orbitals arc included. Re- 
cently, a mechanism of the OSMT based on the lifting 
of orbital degeneracy is proposed.— According to this 
mechanism, in a system with more than two orbitals, if 
the orbital degeneracy is partially lifted by a crystal field 
splitting, an OSMT may take place even when the band- 
widths are equal. It is suggested that an OSMT trig- 
gered by this mechanism may exist in iron pnictidesi^ 
From Fig. [S] we find that for J > 0, the bandwidth of the 
xy orbital gets stronger rcnormalization than the others. 
This effect becomes more pronounced for larger J. How- 
ever, no OSMT is observed up to J /U ~ 1/2. This is not 
too surprising. On the one hand, the above mechanism 
is proposed by assuming zero inter-orbital hopping. The 
four-orbital model discussed in this section contains non- 
zero inter-orbital hoppings, which enhance the orbital 
correlations and favor an one-stage Mott transition. On 
the other hand, according to the mechanism, an OSMT 
is easier to be realized if the orbital whose degeneracy 
is lifted is at half-filling but the degenerate orbitals are 
away from half-filling, so that the lifted orbital is local- 
ized while the degenerate ones are still metallic. In the 
four-orbital model we consider, the degenerate xz and 
yz orbitals are exactly at half filling, making the Mott 
insulator more stable. 



V. CONCLUSION 

In conclusion, we have studied the metal-insulator 
transitions in several multi-orbital Hubbard models for 
the parent iron pnictides using a slave-spin formulation. 

In the two-orbital model, a transition to a Mott insu- 
lator generally exists at half-filling. The Hund's coupling 
reduces the critical coupling significantly. We find that 
the nature of the Mott insulator depends on the ratio 
U' /U. For U' < U, the insulator is a high-spin Mott 
state with zero double occupancy but a spin triplet. For 
U' > U, by contrast, the insulator is an orbital-Mott 
state with spin singlet and finite double occupancy. The 
low-spin orbital-Mott state is unstable to a band insula- 
tor if the orbital degeneracy is lifted, and can be viewed 
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as a special case of an orbitally polarized insulator. 

The phase diagram for the metal-to-insulator transi- 
tion in the more realistic four-orbital model contains ad- 
ditional features. We find a transition to a high-spin 
5 = 2 Mott insulator when the Hund's coupling is strong 
enough {J/U > 0.009). This high-spin Mott insulator 
is an analogy to the S = 1 Mott insulator in the two- 
orbital model. At weak including zero Hund's couplings 
[J/U < 0.009) we find a transition to intermediate-spin 
insulating states. Such a transition is an orbitally selec- 
tive metal-insulator transition, namely, different orbitals 
undergo transitions to different insulating states at a sin- 
gle critical value Uc- At J = the transition leads to a 
band insulator for two (xy and — y^) orbitals and a 
Mott insulator for the other two {xz and yz) orbitals. For 
nonzero J in this regime, the transition in the xz and yz 
orbitals is to the Mott insulator, but an orbitally polar- 
ized insulator is stabilized in the xy and — orbitals 
on the insulating side due to the pair-hopping term. As 
J is increased further, the intermediate-spin orbitally po- 



larized insulator undergoes a low-spin to high-spin tran- 
sition to a high-spin Mott insulator. 

The existence of a Mott transition in the multi-orbital 
models with an even number of electrons per Fe provides 
the theoretical basis for the recent finding of a Mott insu- 
lator in the iron oxychalcogenides with an expanded Fe 
lattico^i. In addition, it strengthens the notion that the 
iron pnictides are located in proximity to a Mott local- 
ization transition. 

Note Added. After this paper was first submitted for 
publication and posted on the arXiv preprint listing, sev- 
eral studies also discussed the relation between the band 
and orbital pictures in related contexts^i*^, and a num- 
ber of works used the slave-spin or a related slave-roter 
method to investigate the metal-insulator transitions in 
related multi-orbital models and systems^"—. 
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